Computing all 5-integral points on elliptic curves 

Attila Petho \ Horst G. Zimmer, Josef Gebel and Emanuel Herrmann* 



1 Introduction 

Let a,b G Z with Aq = 4a^ + 27fo^ ^ 0. Define a finite set of places S = 
{qi, . . . , Qs-i, Qs = cxd} of Q and put Q = max{gi, . . . , Qs-i}- Consider the group 
S(Q) of rational points or Mordell-Weil group of the elliptic curve 

E : y'^ = +ax + b (1) 

over Q with discriminant Aq and absolute invariant j = 12^^. The multi- 
pHcative height of a rational point P = {x, y) G i?(Q) is defined as the following 
product over all places g of Q (including q = oo): 

H(P) = l[ma^{l,\xU}, 

Q 

where the \x\q^s are the normalized multiplicative absolute values of Q corre- 
sponding to the places q and satisfying the product formula. 
Let E{Iis) denote the set of 5— integral points of S(Q), i.e. 

E{Zs) = {P= {x,y) G E{q)\Hs{P) < 1}, 

where 

HsiP) = nmax{l,|z|J. 

Siegel [22] proved in 1929 that the number of integral points on an elliptic curve 
E over an algebraic number field IK is finite, and Mahler [18] generalized this 
result in 1934 to 5-integral points. Moreover, beginning with the pioneering 
work of Baker [1] , several authors derived bounds for the size of the coordinates 
of integer points on elliptic curves E over K. Coates [5] proved that E{Zs) 
contains only finitely many elements and that these can be effectively computed. 
The proofs used algebraic number-theoretical methods; the equation (1) was 
transformed into finitely many S— unit equations. 

A completely different method for proving the finiteness of E{Zs) was proposed 
by S. Lang [16] and Zagior [.31]. Here one makes use of the group structure of 
E{Q) and invokes properties of ordinary and g-adic elliptic logarithms. After S. 
David [7] established an explicit lower bound for linear forms in complex elliptic 
logarithms, Lang's idea could be transformed into an algorithm for computing 
all integer points on elliptic curves. This was done by Gebel, Petho and Zimmer 
[10] and independently by Stroeker and Tzanakis [27]. The algorithm of Gebel, 
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Petho and Zimmer was implemented by Gebel in the computer algebra package 
SIMATH and is available free of charge. 

N. Smart [25] generalized the above method to an algorithm for finding all S- 
integral points on elliptic curves. Unfortunately, his approach, which seems to 
work well in practice, depends on an unproved lower bound for linear forms 
in (/-adic elliptic logarithms. If the rank of i?(Q) is at most two, then such a 
bound was obtained by Remond and Urfels [21] and applied by Gebel, Petho 
and Zimmer [11] to find all S-integral points on Mordell's curves = + k, 
with \k\ < 10^ and such that the rank of the curve is strictly less then three. 
In this note wo combine the advantages of the methods of Siegel-Baker-Coates 
and of Lang-Zagier to overcome the absence of an explicit lower bound for 
linear forms in g-adic elliptic logarithms. This was made possible by a recent, 
completely expHcit upper bound for the S"— integral solutions of (1) established 
by Hajdu and Herendi [15]. 

In Section 2 we state our theoretical result and in the subsequent section we 
prove it. In Section 4 we deal with the computational aspects of the main result 
and in Section 5 we explain how we compute g-adic elliptic logarithms. Finally 
we show through an example how our method works in practice. 

2 Notations and the Theorem 

To state our result we require several definitions. Let Pi, . . . ,Pr denote a basis of 
the Mordell-Weil group and let g be the order of the torsion subgroup E{<Q)tors 
of S(Q). We have < 12 by a famous result of Mazur [20]. Let h denote the 
Neron-Tate height as well as the associated symmetric bilinear form on i?(Q). 
Designate by A the smallest eigenvalue of the positive definite regulator matrix 

ihiPi,Pj))l<iJ<r. 

Let p(m) be the Weierstrafi p— function corresponding to the curve E{C). Let 
be its fundamental lattice and oj its real period. There exists, for any P = 
{x,y) G E{C), an element u € C/il such that {x,y) = {p{u) , p' (uj) . This u 
is called the (complex) elliptic logarithm of P. In the sequel Ui^oo denotes the 
elliptic logarithm of Pi, i = 1, ... , r. We put u[ ,^ = gui^oc/'^- 
For a prime q o{Q, E denotes the reduced curve E modulo q. Let J\fq = #S(Fg ) 
be the number of rational points on E/¥q . With the order g of the torsion group, 
we define the number 

m = nig = lcm{g,Afq). 

Finally, for the finite places q £ S, let u'^ ^ denote the g— adic elliptic logarithm 
of mPi for i = 1, . . . , r. Now we are in a position to state our main result. 

Theorem Assume that the S— integral point P = {x,y) £ E{Iis) has the 
representation 

r 

P = Y,niPi+T (2) 

8=1 
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with n, G Z, i = 1, . . . , r, and T a torsion point of E(Q). For N = max{|nj|,i = 
1, . . . , r}, we have 

N<No = V(ci/2 + C2)/A, (3) 

with 

C2 = logmax{(8|a|)i/4,(32|&|)i/8}, 

Ci = 7 • 1038^+«s20.+ 15Q24(iQg* Q)^^-'c3(l0gC3)2(c3 + 20(s - 1)C3 + log(eC4)), 

where log* Q = maxjlogQ, 1} for Q = max{gi, . . . ,9s-i}, 

32 / I 

C3 = y^A^(^8 + -log|Ao| 

C4 = 10^max{16a^256i/|A^''} 

with Ao = 4a3 + 27b'^. 

Moreover, there exists a place q G S such that 

<C5exp{-{A/s)7V2+c2/s}, (4) 

Q 

with rir+i £ Z i/g = 00 and n^+i = otherwise, and with C5 = y/Sg/to if q = 00 

and C5 = 1 otherwise. 

Remarks 

1. If r < 2, then we can use the estimate of Remond and Urfels [21] for hnear 
forms in g-adic elHptic logarithms to derive an upper bound for N (see 
[11]). The example shown in Section 6 reveals that this bound can be 
much smaller than our bound A^o- 

2. The proof of the Theorem will show that, in the Proposition below, Ci 
can be replaced by any upper bound for the height of 5— integer points in 
i?(Q). For example, in the case of s = 1, one can take the smaller constant 

C'l =5-10«4c3logC3(c3+logC4)) 

instead of speciaHzing to s = 1 in the formula for ci (see [15], Theorem 
1). Then we obtain 

iV < = ^(c;/2 + C2)/A 

instead of (3). 

In the case of s = 1, Gebel, Petho and Zimmer [10] * as well as Stroeker and 
Tzanakis [27] estabhshed another bound for N, namely 

iV < TVi = 2''+^y^\og^''+^^/\c7{r + 3y+^) 

*The constants given here differ from those in [10] because in [10] (on pages 179-180) 
two corrections are necessary: First r must be replaced by r + 1 and second, C requires an 
adjustment to the new version of the Theorem in [7]. 




E 

i=l 
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with 



f log(2^/2^g/c.) 1 
C6 = max^ J ,1V, 

cr = max I p 10^1 ^ Hl^S^'' 

h = logmax{4|aj2|,4|6j2|, |ii|} for j = — , 
h,h e gcd(ji, j2) = 1, 

and 

log V, > max { hiE,), K '"'^b''"" } ' « = 1, ■ ■ ■ , 
logVo > max{/i, 

It seems hard to compare the size of A^q and Ni because they depend on different 
parameters of the curve. The numerical examples of Section 6 suggest that 
TVi < Nq if r < 2, but A^i > Nq if r > 2. Therefore, it makes sense to compute 
both bounds Nq and A''i before starting the reduction procedure, which is based 
on inequality (4) and does not depend on the way, by which an upper bound 
for N was obtained. 



3 Proof of the Theorem 

The proof of the Theorem is divided into three parts: First we obtain an explicit 
upper bound for the ordinary height h{P) of the 5— integral points P in S(Q). In 
the second part we prove a relation between h{P) and TV and derive (3). Finally, 
we prove (4) by using properties of complex and g— adic eUiptic logarithms. 



3.1 Estimation of /i(P). 

Proposition Let E{Q) be the elliptic curve (1) and S he the finite set of rational 
places (containing the infinite place oo) defined above. Then all S— integral 
points 

P=ix^y)= (^-^'^y tvX^"^, gcd(4,C) = l = gcd(r?,C) 
in E{Q) satisfy 

max{h{x),h{y}} < Ci, 
with the constant Ci defined in the Theorem, where 




= logmax{|a|, \f3\}, a,f3 £ Z,gcd(a,/3) = 1. 



Proof: [15], Theorem 2. 
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3.2 Height estimates and the proof of (3). 

Let us start out from an elliptic curve 

E:Y'^+ aiXY + aaY = + a-iX"^ + a^X + ag (a^ £ Z) 

in long Weierstrafi form over Q. Tate's quantities are 

b2 = af + 4a2, &4 = 2a4 + aiOs, bg = al + Aa^, 
bg = a'lae + 4a2a6 — 010304 + 0203 — 04, 
C4 = foi - 2464, ce = -bl + 866264 - 21666- 

E has discriminant 

A = -6^68 - 8bl - 27bl + 9626460 = -I6A0 

and absolute invariant 
3 

The invariant differential is 

dx dy 



2y + aix + 03 + 2a2X + 04 — aiy 

We begin by estimating heights. In the introduction we defined the multiplica- 
tive height of a rational point P = {x,y) g ^^(Q). The ordinary additive height 
is then 

hiP) := i log HiP) 
and the Neron- Tate or canonical additive height is 

HP) := hm 

From Theorem 5 in [34] (cf. also [24], [32], [33]), one readily derives the inequal- 
ity 

HP) > HP) - 4 for p G E(q) 

with the constant ^ 

C2 := 2 + fi^), 

where ^ ^ ^ 

/ioo := log max{|62|oo, \b4\S0, Ibelio, Ibs]^}- 

We now specialize this height estimation to the curve in short Weierstrass form 
(1). We have oi = 02 = 03 = and 04 = a,ae = 6, so that yuoo becomes 

/ioo = logmax{(2|a|)i/2, {4\b\)'/\ \a\'/'} = logmax{(2|a|)i/2, {4\b\)'/'}. 

Hence 

HP) > HP) - C2 for P e E(Q). (5) 
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with the constant ci defined in the Theorem. 

Let, as before, Pi, . . . ,Pr be a basis of the free part of the Mordell-Weil group 
E{Q). Then P e E{Q) has the (unique) representation (2) with rii e 1^,1 = 
1, . . . r, and T a torsion point. Put TV = max{|ni|, ? = 1, . . . , r}. The canonical 
height h can be extended to the r— dimensional real vector space i?(Q) ®z K 
and the extension of /i is a positive definite quadratic form on this space. 
As explained in [10], this leads to the lower estimate on non-torsion points 

h{P) > XN"^ 

involving the smallest eigenvalue A of the regulator matrix associated with the 
basis points Pi, . . . ,Pr £ -E(Q) via the canonical height h. This, together with 
(5), implies that 

hiP) > - C2 for P e EiQ) \ EiQ)tors ■ (6) 
Let now P = {x,y) £ E{Zs)- Then, since h{P) = ^\ogh{x), we obtain 

HP) < Ci/2 

by the Proposition, with Ci as defined in the Theorem. On combining the lower 
and upper estimates for h{P), we obtain 

-C2 < ci/2. 

This impHes the asserted inequality (3): 

iV< V(ci/2 + C2)/A) 

3.3 Elliptic logarithms and the proof of inequality (4). 

Let now P = {x, y) £ be an S'-integral point and choose q & S such that 

\x\q =max{|a;|,i,...,|a;|,._i,|a;|oo}. (7) 
In the sequel we keep this q fixed. First we conclude that 
H{P) < \x\l for s := 

and hence that 

HP) < I log \xU. (8) 
Combining (6) and (8) yields (cf. [8]) 

< cs exp {-CbN^) (9) 

\xW 

with 

fC2\ A 

C8 := exp I — , Cg := -. 

\ s / s 
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To transform the upper bound for \x\q '■' into the upper bound for the left hand 
side of the inequality (4), we must study the properties of elliptic logarithms. 
Two cases are to be distinguished, the complex case of g = oo and the g-adic 
case. 

Case 1: q = oc £ S. 



Let p(m) be the Weierstrafi function associated to the curve (1). Let fl = {lo, L02) 
be the period lattice of p{u), generated by the two fundamental periods w and 
0J2, where w > is real and 0J2 is complex. Let P = {x,y) G i?(Q) and denote 
hy u = u{P) its elliptic logarithm, i.e. that residue class u G C mod fl for which 
P=ip{u),p'iu)). 



Then we have 



(see [31]). 
Let 



f 



(mod 0) 



Vt^ + at + b 



+ at + b={t-a){t- /3){t - 7) 
be the decomposition of the cubic polynomial into linear factors with roots 

a < /3 < 7 if a,/3,7 e M 

and 

a = P eC, J eW otherwise. 

Define 



M 

and put 



, if 7 > 

exp /ioo . . ^ „ 
^TTJ^I'if 7<0 



xq = max{a + /3, 2j} + M. 
With this notation, we have the 

Lemma Let a; £ IK such that x > max{0,a;o}. Then 

dt 

< 



f 



Vt^ +at + b Vx' 

Proof: [10], pp. 176 and 181-183. Note that in [10], the number x is assumed 
to be the first coordinate of an integer point of i?(Q), but in the proof this 
assumption was not needed. 

Let now P = {x,y) be an 5— integral point on S(Q) with first coordinate of 
absolute value 

|a:|oo > max{a;o,expjUoo}- (10) 

Denote by u its ehiptic logarithm. Assume furthermore that (7) holds with 
q = 00. Then, by the Lemma and by (9), we obtain, on normalizing u to 

< |w|oo < 



<V&C8 expi-CgN"^). (11) 
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Let Ui^oo = u{Pi),i = 1, . . . , r, denote the elliptic logarithms of Pi, ... , Pr- Then 

(r \ r 

i=l J i=l 

where g is the order of the torsion subgroup of S(Q). Moreover, there exists an 
integer n^+i £ Z such that 

|"i(5Wi,oo) + • • • + nr{gUr^oo) + nr+iU}\oo < Vs g Cg exp{-CgN'^). 

Dividing by lo we obtain (4) in the case of g = oo. 

Case 2: q = e S (for some i £ N such that 1 < i < s — 1). 
Here we use g-adic elliptic logarithms (cf. [8], [23], [25]). 

For the convenience of the reader, we explain in some detail how one proceeds 
in the g-adic case. Let be the g-adic completion of Q and be its ring of 
g-adic integers. Denote by 

EiiQg) ■.= {PGEiQg) \ P = d} 

the kernel of the reduction map modulo q, where E is regarded as a curve over 
Qg and P, O are the reduced points P,0 modulo q. Designate by £{q'Lq) the 
formal group associated to E (cf. [23], [25]). We consider the isomorphism 

f(gZ,) Ei(Q,) 

J if;z = 1 



where 



--, w{z) -- 

y y 



The equation for w = ■w{z) inferred from the long Weierstrafi equation for -E(Q) 
becomes 

w = z^ + {aiz + a2z'^)w + (03 + a4z)w'^ + aeiu'^ = f{z,w). 

A recursive procedure based on this equation (see [23]) leads to the power series 

w = z^ + aiz'^ + (af + a2)z^ + (af + 2aia2 + a^)z^ 
+ {a\ + 3ala2 + SaiUa + a| + a4)z'^ + ■■■ 
e Z[ai,a2,a3,a4,a6][[z]\. 

This is the unique power series in z satisfying the relation 

w{z) = f{z,w{z)). 

From it we also get the Laurent series for x and y, viz. 

^i^) = ^ = - ^ -a2- a-iZ- {a4+aia^)z'^ , 
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The invariant differential has the expansion 

Lo{z) = (1 + aiz + (af + a2)z'^ + (af + 2aia2 + 03)2;^ 
+ {af + 3afa2 + 60103 + al + 204)2;'' + ■ • ■)dz. 

Note that in these expansions the coefficients of the powers of z each have the 
same weight depending on the exponent of z. 

The g-adic elliptic logarithm is now the image under the homomorphism to 
the additive group Ga (over the completion of the algebraic closure of Qg.) 
defined as follows: 

P = ix,y) ^ ■ipg{P)=Jujiz) = z+^z'^ + fz^ + ---. 
In particular, the g-adic logarithm has the properties 

and 



As before, let E be the reduced curve E modulo q, and denote by Afg = j|S(Fg) 
the number of rational points on E/¥g. With the number 

m = lcm((7, Ag) 

defined in Section 2, we have from the Lutz filtration of E (see [17]): 

mPi=:P/e£i(Q,) {i = l,...,r) 
for the generating points Pi of E{Q), and 

mPr+i = O 
for the torsion points Pr+i £ E{^tors- 

The representation (2) of an 5-integral point P = {x, y) € E{Q) then gives rise 
to the representation 

P' =n[Pi+--- + nlPr = niPl + ---nrPl. {n'^ := mrii G Z) (13) 

of its m-multiple P' = {x',y') = mP € -Bi(Qg). In analogy to (12), we have the 
Laurent series 

I 1 Oi I , . ,2 

— ~ ~2 ~ ~ ~ "2 - aaz - (04 + 0103)2 - ■ ■ ■ 



'w{z') z' z 
and this expansion entails the estimate 
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where we use the abbreviating notation u'^ := tjjq{P') for the eUiptic logarithm 
of P'. 

On combining the inequalities (9) and (14) and observing that \x'\q > \x\q, we 
end up with the upper estimate for the g-adic elhptic logarithm u'g = il^q(P') of 
the point P' = {x',y') = mP = m{x,y): 

\u'qU < — |t7^ < -jrj^ < cs exp(-C9A'2). 

\X \g \X\g 

From (13), we have the relation 

= n[ui^g + h n'^Ur,q = niu[^^ + . . . + riru'^^^ {n[ = mrii e Z) 

between the elhptic logarithms u'q = tjjq{P') of the m-multiple P' = mP of 
the given 5- integral point P, Ui,q = ipq{Pi) of the generating points Pi and 
u!i g = ipgiPj.) of their m-multiples P/ = mPj g E{Q), and inequality (4) is 
proved in this case too. 

4 Computational aspects of the Theorem. 

The Theorem, together with some computational techniques for elliptic curves 
and with diophantine approximation methods, in practice makes it possible to 
compute ah S"— integral solutions of (1). To this end, we must carry out the 
following steps: 

(i) Compute a basis P\,. . . ,Pr of the Mordell-Weil group E{^. This can be 
achieved, for example, by employing the algorithm [9] (see also [19], [33]). 
Having obtained the basis, one can easily compute the constants Ci , C2 and 
A and thus the upper bound A^q for N , defined by the inequality (3). 

(ii) Nq is usually very large (see the examples), but we also have the estimate 
(4) at one of the places q G S. Some observations on the computation 
of rational approximations of the elhptic logarithms with the necessary 
very high precision will be made in the next section. The solutions of the 
system of inequalities (3) and (4) are then found by the reduction method 
of de Weger [29]. Since this method is well-described in [25], we omit the 
details here. Of course, one must perform de Weger's reduction for all 
places q e S. Let Mq denote the bound for N after the reduction with 
respect to a prime q £ S, assuming that (7) holds for this q £ S. Then 
we change the value of Nq to ina.x{Mq,q G S}. If possible we iterate the 
reduction. 

(iii) Experience shows that the last value of A^q obtained after several reductions 
is sufficiently small in order to enable one to test the points Yll=i niPi + T 
within the range TV = maxjjnjj, « = 1, . . . , r} < A^o and with T g EtorsiQ) 
for S'-integrality. 
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However, this can be a lengthy procedure. Therefore, we describe here 
some ideas as to how to proceed in a more economical way. If the r— tuple 
n = (rii, . . . , rir) of integer vectors is such that TV = max{|nj|, i = 1, . . . , r} < 
A^o, then we first test for q G S hy fioating point approximation (if q = oo) 
or by g— adic approximation (if g is a finite place) of whether or not 
(4) is true. If, for an n, inequality (4) does not hold for some q £ S, 
then there is no S"— integral point on the curve given by the vector n in 
the representation (2). Otherwise we must test all linear combinations 
Yll=i n-iPi + T with varying torsion points T for 5— integrality. 

(iv) There is one extra search to be carried out eventually. This is due to 
the fact that, if (7) holds for q = oo, then (11) and (4) hold only for 
those 5— integral points P = {x,y), whose first coordinate satisfies the 
inequality (10). We must also determine the S'-integral points for which 
(10) is violated. Then, by (7), we have 

\x\g. < max{a;o,exp/ioo}, i = l,...,s-l. 

Write 

x=—, — ^ a.-i with gcd(^,gi) = 1, i = l,...,s-l. (15) 
This yields 
so that 

a, < logmax{a;o,exp/ioo}/loggj, i = l,...,s-l. (16) 

Hence the exponents are bounded. We note that the exponents 
are even. The negation of inequality (10), together with the boundedness 
relations (16) for the coordinates (15), shows that there are only finitely 
many 5-integral points, which are to be tested in this case. 

5 Computation of the g-adic elliptic logarithms 

The complex eUiptic logarithms can be computed either by applying the arith- 
metic-geometric mean method (see H.Cohen [6]) or by employing a method of 
Zagier [31]. 

In contrast to the case of complex elliptic logarithms, there exists, for the com- 
putation of g-adic elliptic logarithms, a quadratically convergent procedure only 
for Tate curves (see [14]) *. That is why we must apply a different g-adic pro- 
cedure here. In fact, we simply use the formulae of Section 3. However, we 
must overcome some technical difficulties if we need to calculate the g-adic el- 
liptic logarithms to a precision of up to 1000 or more digits. Of course, we can 
compute the formal Laurent series expansion of w(T) by Newton iteration. How- 
ever, unfortunately the coefficients in this expansion soon become very large. 

*We wish to thank J.H.Silverman for calling our attention to the paper [14] 



11 



We circumvent this problem by invoking the Lutz filtration of E over and 
calculating g-adic approximations to a sufficient high precision as follows. 
Suppose that, for some given positive integer t, we have computed an approxi- 
mation 



^(T) (mod T*+i) 



of the formal elliptic integral ip{T) = J w{T) dT, where di = 1. We want 

to compute an approximation tpg{P) of the g-adic elliptic logarithm tpg{P) of 
P 6 -Ei(Q,) such that 

IMP) - UP)U < q-" 

for some positive integer n. (In practice, t is fixed and n depends on q and 
is much larger than t.) Put V = [j] and define Py = q^P- Then we have 
Pv G Ev+iiQg) by the Lutz filtration theorem [17]. Now we can compute a 
suitable g-adic approximation Py of the point Py since E is considered as a curve 
over <Qg as in Section 3.3. Here a suitable approximation means that, for the 
point Pv = {xp^,ypy) e EiQg), we choose Py = (xp^ /q''^^+^\yPv /q^^^+^^) 
with xpy,ypy e Z such that 



XPy 

ypv 



^Py 

ypv 



< 9"" 



with n' = n + V. To this end we take zpy G Z such that 

\zPv\g 
Then we have 



^Py qV+l 


and 


, XPy 

zpy + - — q^ 


< q~" 


ypv 


Q 


ypy 


Q 



ZPy + 



XPy 

ypy 



<q-\ 



which implies 

i'izpv) - i' 



XPy 

ypy 



On the other hand, we have 



XPy 

ypy 



< 



XPy 

ypy 




< 



H+1 



t + 1 



XPy 

ypy 



t+1 



i 

i=t+l 

-(i+l)(v+l) 



XPy 

ypy 



< 



t + 1 



< 9"" 
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by the choice of V. Thus 



ypv 



Q 



We have 



hence 



ypv 



\MP)-^iSpv)/Q''U<q-''^=q 



i.e. 4'{Spv)/q is the required approximation of -0g(P), the g-adic elliptic loga- 
rithm of P. 

6 Comparison of the Theorem with earlier esti- 
mates 

We first consider the S = {3, 5, oo}— integer points of the curve 

EsmiQ) -.y^ + y = x'' - 2x, 

treated by Smart [25]. He proved that Pi = (0,0), Pj = (liO) form a basis of 
the Mordell-Weil group of Esm over Q and that Esm has no non-trivial torsion 
points. Moreover, he showed that if 

P = niPi+n2P2 

is an S— integral point on Esm with max{|ni|, \n2\} < W^^, then ±(711,712) = 
(0, 1), (0, 2), (0, 3), (1, 0), (1, 1), (1, -1), (1, -4), (2, 0), (2, -1), (2, 2), (2, -2), 
(2, -3), (3,0), (3,1), (3, -2), (4,-1), (1,2). 

Using the transformation X = 36x -I- 12,y = 216y + 108, we obtain the short 
Weierstrafi form 

E'smiQ) :Y'^ =X^ - 3024X -I- 46224. 

It is clear that if P = (a;, y) is an 5— integral point on Esm, then P' = {x' ,y') = 
{36x + 12, 216y + 108) is an 5— integral point on E'g^. Thus it suffices to find 
all 5— integral points on E'g^. Let P' = {x',y') be an 5— integral point on Eg^ 
with 

P' =niP'i +n2P'2 
for P'l = (12, 108), P'2 = (48, 108). Then the Theorem implies that 

TV = max{|ni|, \n2\} < Nq = 6 ■ 10^^\ 

as follows by a simple computation. This bound is much larger then the bound 
used by Smart. However, by performing one de Weger reduction with this 
unconditional bound, we obtain 

max{|ni|, \n2\} < 106. 
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This shows that Smart indeed found all 5— integral points on Esm- 

As the rank of Esm is two, one can also apply the theorem of Remond and 
Urfels [21] to estimate N in the way described in [11]. Using that method one 
easily verifies that 

TV < 4.7- 10"^ 

Here the worst case appears when maxljitja, |m|b,|m|cx3} = jwjoo- This is a much 
smaller bound for N than the one which follows from the Theorem, and hence 
we obtain another proof that Smart found all 5— integral points on Egm, this 
time without any further computation. 

In Remark (2) and thereafter, we gave two different bounds TVq and Ni for N in 
the case s = 1. Now we are going to compare them for some curves, for which 
all integer points were found by the eUiptic logarithm method. 
The first set of examples is taken from [13], where the classical Mordell equations 

y"^ = +k 

were solved t for jfcj < 10^ In Table 1 wo list, for some fc's, the rank rk of the 
curve y"^ = + k and the values of A^i and Aq. 



Table 1 



k 




iVi 




108 


1 


2.8 • 102« 


1.8-10" 


225 


2 


1.3- 10"! 


4.5 • 10" 


1025 


3 


5.5 ■ 10*5° 


3.1 ■ lO''^ 


2089 


4 


1.1 ■ 10*^ 


7.7-10^2 


-28279 


5 


2.1-10112 


2.0 - 10^2 



Next, in Table 2, we take some examples from [3], where the family of eUiptic 
curves 

2y2 = a;3 _ 36a. _ 864fc(fc - l)(2fc - 1) 

was considered. 



Table 2 



k 


Tk 


iVi 




% 


1 


1 


8.9 


1023 


1.3 


10^1 


3 


2 


5.8 


1039 


1.8 


1044 


7 


3 


2.4 


1060 


5.9 


10^5 


20 


4 


2.1 


1086 


2.9 


10^7 



tThe rank one cases left open in [13] have been settled by K. Wildanger in his PhD thesis 
"Uber das Losen von Einheiten- und Indexformgleichungen in algebraischen Zahlkorpern mit 
einer Anwendung auf die Bestimmung aller ganzen Punkte einer Mordellschen Kurve" , Berlin 
1997. In those cases, the Mordell curves have no integral points - as conjectured by the 
authors. 
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Finally, we consider four curves from [8], [10] and from [28] of the form 

y'^ = + ax + b. 

In Table 3 we exhibit, for certain coefficients a, b, the rank r of the curve defined 
by (1) and the corresponding values of A'^i and A^q- 



Table 3 



a 


h 


r 


N, 




-203472 


18487440 


5 


2.3 • 10"^ 


7.9 - 10*'' 


-1642032 


628747920 


6 


1.1 • 10"* 


2.1 - 10*9 


-147952 


21183760 


7 


2.7-101" 


1.1-10*7 


-5818216808130 


5401285759982786436 


8 


8.67-10224 


2.33 - 10^^ 



The examples indicate, that A'l < TV,^, if r < 2, but Ni > Nly, if r > 2. Of 
course, for rank curves it is quite easy to find all integer points by examining 
the algebra of the curve. Our comparison shows that Ni is significantly smaller 
than iVg if r = 1. For r = 2, Ni and Nq are of a similar size, but Ni appears to 
be slightly smaller than Nq. For r > 3, the constant A^q becomes much smaller 
than A^i, and with growing rank r, the difference becomes rather significant. 
This can be very well observed in the above rank 7 and 8 examples. The reason 
for this phenomenon probably lies in the fact that in David's [7] lower bound, 
which enters into the bound A^i, the constant 

C = 2.9 - 10^(''+2)42('-+i)'(r + 2)2''^+i3'-+23.3 

shows up, and this constant depends exponentially on r^. On the other hand, 

the constant term of A^'q docs not depend on r. 

Our examples make it plausible that, in order to find ah integer points on elhptic 
curves, it makes sense to first compute both bounds A^i and A^q and then start 
the reduction procedure with the smaher one. 

7 An example 

To illustrate our method, we choosen the rank-four curve 

E: y'^ =x^ - 172x + 505 (17) 

from the paper A. Wiman [30] (see also Bremner and Tzanakis [4]). Wiman 
showed that there are at least 58 integer points on E. Here we prove that there 
are exactly 58 integer points and, moreover, 144 S- integer points on E, where 
S = {3, 5, 7, oo}. More precisely, we prove the following theorem. 
Theorem' Let S = {3,5,7, oo}. Then all S -integer points P = {x,±y) with 
X = ^,y = C S Z,gcd(^, C) = gcd(7;,C) = 1 on E over Q are those 

displayed in Table 5. 

Proof. First we compute the quantities 01,02,03 and 04 of the Theorem in 
Section 2 to obtain A^o in dependence on A only. For this purpose, we do not 



15 



need any geometrical information about the elliptic curve E. We have a = —172 
and b = 505, hence 

Ao = 2198992, 02 = 1.81, C3 = 8.7 • 10^ and a = 8.35 • 10^^ 

Taking s = 4 and Q = 7, we obtain 

Cl = 4.564- 10^°^ 

and this yields 

A^o = 4.777 ■ 10l^V^/A. 

Now we consider E as an elliptic curve over Q. Inserting the coefficients of E 
into SIMATH we obtain, after some minutes of computation, that E has trivial 
torsion, its rank is 4, and a basis of its Mordell-Weyl group is given by the 
points: 

Pi = (12,13), P2 = (-14,13), P3 = (-1,26) and ^4 = (38,221). 

The regulator of the curve is R = 2.79532, the real period is w = 0.808974 and 
the least eigenvalue corresponding to the chosen basis is A = 0.7467531. 
We have g = 1, hence = Mq, and the actual values of are ms = 7, 7715 = 10 
and my = 12. By the Theorem, for each q G S, we have to solve the following 
diophantine approximation problems 

< exp{-0.149- +C10}, 

Q 

TV = max{|ni|,...,|n4|} < iVo = 5.53 • lO^^^ 

for each q & S, where cio = 1.62, if g = oc and cio = 0.362 otherwise. 
In order to reduce the huge upper bound for N, we first take q = 00 and 
perform de Weger reduction with C = 10^^°. We obtain the new upper bound 
N' < Moo = 69. Next, for each q G S \ {00}, we compute the g-adic elliptic 
logarithms oi Pi,i = 1, . . . , 4, with precision at least 

ris = 1281, rir, = 875, riy = 723. 

This precision is necessary in order to carry out the g-adic de Weger reduction. 
For this purpose, we use the method of Section 5. We have chosen t = 200 and, 
for safety reasons, 

1/3 = 8, ^5 = 6, V7 = 5, 

and computed the approximations of g^'Pj for q € {3,5,7} and i = 1,...,4 
with 

na = 1700, rig = 1300, n'j = 1100 

g-adic digit numbers. In Table 4 we have listed the first and last three significant 
digits of tpqlniq Pi), i.e., the digits 02, 03, 04, an^_2, an,_i and a„^ assuming, that 

i>q{mqPi) K ^aiq\ 

Note that, as m^Pj € i5i(Q), we always have Oi = 0. 
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Table 4 



i 


1 


2 


3 


4 


Pi 


(12,13) 


(-14,13) 


(-1,26) 


(38,221) 


M^Pi) 


(1,2,1,. ..,0,2,2) 


(2,0,1,. ..,2, 1,2) 


(1,1,0,. ..,2,2,2) 


(0,2,2,. ..,0,1,1) 


MWPi) 


(3,0,4,. ..,3,3,4) 


(4,4,0,. ..,0,3,1) 


(2,2,3,. ..,3,3,0) 


(4,3,4,. ..,0,1,1) 




(0,6,0,. ..,2, 3,4) 


(5, 2,5,. ..,0,2,4) 


(2,1,1,. ..,5,6,2) 


(1,4,4,. ..,1,2,1) 



Now we perform the g-adic de Weger reduction with the values C3 = 3^^^^, Cg = 
5*''^ and C7 = 7''^^ and obtain the new bound 

No = 124 = max{7Woo,^3,^5, A^t}- 

This new upper bound for N can be further reduced. On repeating this reduc- 
tion process 3-times, we eventually get N < 17, which cannot be reduced any 
further. In order to compute all ^-integral solutions of (17) with this last bound, 
we use the sieving procedure described in [8]. Finally, we must also determine 
the S'-integral points for which (10) is violated. In the case under consideration, 
we have 

la;],, < max{22.56, 31.86} = 31.86, i = 1, . . . , 4, (18) 

thus Q!2 < 4 and (23, as < 2, where we used the notation of Section 4, part (iv). 
A simple computation shows that all the solutions of (17) satisfying the above 
inequality (18) are among those already found. Theorem' is proved. □ 

Table 5 



S'-integral points P = {x, y) = 




E: = a;3 - 172a; + 505 
for 5 = {3, 5, 7, 00} 



rank 


4 


basis 


Pi = (12,13), P2 = (-14,13), P3 = (-1,26), P4 = (38,221) 


torsion 


trivial 



# 




V 


c 


F 


(ni,n2,n3,n4) 


1 


-14 


13 






(0,1,0,0) 


2 


-12 


29 






(0,0,1,1) 


3 


-10 


35 






(-1,0,-1,0) 


4 


-8 


37 






(0,-1,0,-1) 


5 


-2 


29 






(1,0,-1,0) 
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E : y"^ = - 172x + 505 (continued) 



41 

it 


c 
4 




C 


TP 
r 


(ni, 712, na, 714) 


b 


— 1 


2o 






(0,0,1,0) 


I 


o 
2 


1 

lo 






(-1, -1,0,0) 


o 

8 


o 



4 






1 1 n n\ 

(-1, 1,0,0) 


n 

y 


1 o 
12 


1 

lo 






/inn ^^ 
(1,0,0,0) 


iu 


14 


on 
29 






/ inn 1 ^ 

(-1,0,0, -1) 


li 


lo 


A 

4o 






/"n 1 1 f\\ 

(0, -1, -1,0) 


1 o 


1 n 
19 


d4 






(1, 1, 1, 1) 


1 o 


24 


1 ni 
lUl 






/n 1 1 r\\ 

(0,1,-1,0) 


1 /I 
14 


O 1 

ol 


1 CO 

loo 






/inn 1 ^ 

(1,0,0,-1) 


15 


o o 
oo 


oo 1 
221 






/n n n 1 ^ 

(0,0,0, 1) 


io 


YD 


000 






/ n n n^ 

(-2,0,0,0) 


17 


nn 

90 


OAT 

845 






/n 1 1 1 \ 

(0,-1,-1,-1) 


18 


131 


"1 /I no 

1492 






/I 11 1 ^ 

(-1,-1,1,-1) 


19 


lo8 


O "1 'TA 

2171 






/ill n\ 

(1,1,1,0) 


on 


O O /I 

2o4 


/I 'TO 1 

47ol 






/n n n\ 

(0,-2,0,0) 


zi 


/lie 
415 


o /( cn 






/o n n 1 \ 

(2,0,0, 1) 


oo 

22 


4d7 


"1 nno o 

10088 






/n 1 1 1 \ 

(0,1,-1,-1) 


id 


1U4d 


o o OO^ 

66S27 






/ 1 n 1 \ 

(1,0,2, 1) 


O/l 


1200 


A Kr\A 

4oU4o 






(-1, 1, 1, -1) 




lol4 


47o29 






(-1,1,-1,1) 


io 


ono o 
oU28 


1ddo21 






/in n\ 

(1,0, -2,0) 


27 


n cno 

9502 


no£f oo 

926237 






/ on 1 \ 

(-2,-2,0,-1) 


28 


"1 "1 OO o 

1128o 


A A no A f\£i 

1198496 






/ 1 1 1 \ 

(-2,1,1,1) 


on 

29 


1 /inn AC A 

14U24o4 


1 ^f^fnc^T/Jon 

1660877429 






/n 1 o\ 

(0,-3,-1,-2) 


on 


1 oo 
— 12o 


233 


o 



o 



/ 1 n 1 ^ 

(-2,-1,0,1-) 


ol 


4 


559 


o 



o 



/in 1 1 \ 

(-1,0,-1,-1) 


oo 


oo 

28 


1 


3 


3 


/n n 1 1 \ 

(0,0,1,-1) 


33 


106 


287 


3 


3 


/ 1 1 1 1 \ 

(-1,1,1,1) 


O A 

34 


160 


1495 


3 


3 


(-1, -1, 1,0) 


O CT 

35 


oo o 

823 


O O OT /I 

23374 


3 


3 


/n on 1 \ 

(0,2,0, 1) 


00 


1 1 1 O 

11134 


1174769 


3 


3 


/I 1 ^ 
(-1, -2, -2, -1) 


J 






q 


q 



/n n 9 O'i 


38 


-1085 


14680 


9 


32 


(-1,1,2,0) 


39 


-926 


22789 


9 


32 


(2,-1,0,0) 


40 


-536 


26819 


9 


32 


(-1,-1,0,1) 


41 


4612 


305227 


9 


32 


(-1,0,-2,0) 


42 


2182 


811433 


27 


33 


(-1,-2,-1,-2) 


43 


14119 


13113463 


27 


33 


(-1,2,0,-1) 


44 


1709021 


62200333 


81 


3^ 


(2,0,2,1) 


45 


-181 


4628 


5 


5 


(2,1,0,0) 
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E : y"^ = - 172x + 505 (continued) 



# 






C 


r 


(711,712,713,714) 


4o 


TOO A 


o914o7 


r 




r 




/ 1 1 o^ 
(-2, -1, -1, -2) 


47 


1064 


66167 


r 

5 


5 


/O "1 "1 (W 

(2,-1,-1,0) 


A O 

48 


or\ /I 

294 


IdOd 


5 


rr 

5 


"1 on "1 ^ 

(-1,-2,0,-1) 


A n 

49 


3314 


189003 


5 


rr 

5 


/I 1 1 1 \ 

(-1, -1, -1, 1) 


ten 

oU 


OO /I 

— oz4 


one: 1 
z9ol 


r 




rr 



/in 1 1 \ 

(1,0, -1, 1) 


oi 


1 1 m 
1191 


Odd A 

39014 


5 


rr 



/ 1 n 1 ^ 

(-1,0,2, 1) 


52 


34d 


3481 


5 


rr 

5 


/1 1 1 1 ^ 

(1,-1,1,1) 


5d 


— d4 


3809 


5 


rr 

5 


(-1, 1,0, 1) 


54 


— 70>34 


497861 


O CT 

25 


5 


(0,0, -1, -2) 


55 


1754 


1 C\ Al "IT 

1U4117 


O CT 

25 


5 


(0,1,-2,0) 


00 


74Uoo99 


1 CCCO£?01000 

1555zDzl3oz 


^;o c 

o2o 


5 


/1 1 ^ 
(1, -2, -2, 1) 


57 


O O 

i}34 


53677 


15 


3x5 


(0,2,1,0) 


58 


A A'\ C\ AO 

441U4o 


292380031 


A IT 

45 


3x5 


(3,0,0,-1) 


59 


o/? 

—2d 


12571 


7 


7 


/ OH"! (W 

(-2,0,1,0) 


dO 


-657 


coca 

d8dd 


7 


7 


/ 1 01 n\ 

(-1, -2, -1,0) 




16 


7257 


7 


7 


(1,-1,0,1) 


d2 


"1 O O 


or\ A "1 

2941 


7 


7 


(0, 1, 2, 1) 


fad 


cr cr /I 

554 


811 


7 


7 


(1,1,-1,-1) 






91 07 


7 


7 


lu, i, -i, ij 


65 


1082 


29653 


7 


7 


(2,-1,-1,-1) 


66 


140260 


52528645 


7 


7 


(1,0,0,-2) 


67 


-33666 


1488019 


49 


72 


(0,-2,-1,-2) 


68 


-24097 


4109846 


49 


72 


(-1,0,1,-2) 


69 


399988 


252199571 


49 


72 


(1,-2,-2,0) 


70 


-2498 


333593 


21 


3x7 


(-1,2,-1,0) 


71 


28596 


4114769 


35 


5x7 


(2,0,2,0) 


72 


430259 


157313848 


175 


52 X 7 


(-1,3,1,1) 
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